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ABSTRACT 



We present the ensemble properties of 31 comets (27 resolved and 4 unre- 
solved) observed by the Sloan Digital Sky Survey (SDSS). This sample of 
comets represents about 1 comet per 10 million SDSS photometric objects. 
Five-band (u, g,r,i, z) photometry is used to determine the comets' colors, 
sizes, surface brightness profiles, and rates of dust production in terms of 
the Afp formalism. We find that the cumulative luminosity function for the 
Jupiter Family Comets in our sample is well fit by a power law of the form 
N(< H) oc io(°- 49±0 05 ) H for H < 18, with evidence of a much shallower fit 
N(< H) oc io(°- 19±0 - 03 ) H for the faint (14.5 < H < 18) comets. The resolved 
comets show an extremely narrow distribution of colors (0.57 ± 0.05 in g — r 
for example) , which are statistically indistinguishable from that of the Jupiter 
Trojans. Further, there is no evidence of correlation between color and physi- 
cal, dynamical, or observational parameters for the observed comets. 



Keywords: COMETS; PHOTOMETRY; COMETS, COMA 
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1 Introduction 



The physical properties of small planetary bodies offer insight into the forma- 
tion and evolution of the Solar System. However, in situ observational studies 
of the remote regions of the Solar System are limited to the largest bodies 
in the Kuiper Belt, and are still not possible in the case of the Oort Cloud. 
Nevertheless, through delivery of scattered members of their populations, such 
as comets, these regions may be probed observationally. Comets may be dis- 
tinguished from other populations of small bodies due to their activity: the 
production of a gas and dust comae, typically when the comet is at small 
heliocentric distances, due to the sublimation of volatiles. The comet popu- 
lations that are thought to directly sample these two remote regions are the 
Jupiter Family Comets (JFC), from the Kuiper Belt, characterized by gener- 
ally low-inclination prograde orbits, and Long Period Comets (LPC) thought 
to originate in the Oort Cloud, with a large range of orbital inclinati ons, di 



rectio ns, and whose aphelia lie far beyond the orbits of the planets ( IJewitt 



20041 ) . Interactions with the giant planets can give rise to comets that exist 



in-between these two populations, comets with orbits like that of lP/Halley. 
Dynamically, the Tisserand parameter with respect to Jupiter is often em- 
ployed to di stinguish between t he va rious populations of comets based on 
their orbits (jLevison fc Duncanl . 119971 ). Recently a new population of co mets 
has been discovered within the main asteroid belt ( jHsieh fc Jewittl . 120061 ); the 
orbits of these objects are indistinguishable from those of main belt asteroids, 
but they have observed activity in the form of comae and dust trails. 



Building upon the methodolo gy of detecting activ e comets in the SDSS de- 
scribed in our previous work ( ISolontoi et ah] , |2010j ), we present here the anal- 
ysis of properties derived from 35 sets of u,g,r,i, z band phot o metry of 31 



comets observed by the Sloan Digital Sky Survey ( lYork et al.l. |2000|). This 



i nclud es comets found through the methodology described by ISolontoi et al. 
( 120101 ) . along with those recovered through orbital algorithms developed for 
matc hing asteroids found in the SPSS to predicted positions of known aster- 



oids (llvezic et al.l . 120021 ; lJuric et all . 120021 ) 



In Section 2 we present a brief overview of the methods employed to find 
comets in the SDSS. The data and analysis of the observations of both the re- 
solved comets (showing comae), and those that were unresolved (point-source 
like) are discussed in Section 3. We discuss the color, size and dust proper- 
ties of the individual comets, as well as the cumulative luminosity and size 
distributions of the Jupiter Family Comets, in Section 4. 
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2 Observations 



2. 1 Brief Overview of the SDSS 



The data presen ted here are based on th e Seventh SDSS Public Data Release, 



hereafter DR7 (jAbazajian et al.l . 120091 ). which ran through July 2008, and 



contains over 357 million unique photometric objects. Detailed information 
about this data release including sky coverage, changes from previous data 
release s and data quality stati stics can be found at htt p:/ /www.sdss.org/DR7] 
and in lAbazajian et al.l ( 120091 ). Of particular interest to Solar System studies, 
the survey covers the sky at and near the ecliptic from approximately A = 100° 
to A = 225°. The repeat scans of the Southern Galactic hemisphere (Stripe 
82; crossing A = 0°) also pass through the Ecliptic. 



The SDSS is a digital photometric and spectroscopic survey that covered about 
one quarter of the Celestial Sphere in the North Galactic cap and a smaller (~ 
300 deg 2 ) but much deeper survey in the So uthern Galactic hemisphere and be- 
gan standard operations in April 2000 ( see lYork et al.l. boool: IStoughton et al. 



2002 



Abazaiian et^ . l2003l . [2^oll2005[[2009llAdelman-McCarthv et all 12006 



20081 ). SDSS used a dedicated 2.5m telescope ( IGunn et al.l . 120061) to provide 



homogeneous and deep (r < 22.5) pho t ometry in five band passes (Fukug ita et al 



1996 1: IGunn et all 119981 : iHogg et all . 120011 : iSmith et al.l . 120021 : iTucker et al 



20061 ) repeatable to 0.02 mag (root-me an-square scatter, hereafter rms, for 



sources not limited by photon statist ics. Ilvezic et al.l. 120031) and with a zero- 
point uncertainty of ~0. 02-0. 03 mag ( Ilvezic et al.l . 120041 ) . The flux densities of 
detected objects were measured almost simultaneously in five bands (u, g, r, 
i, and z ) with effective w avelengths of 3540 A, 4760 A, 6280 A, 7690 A, and 
9250 A ( Doi et all |2010| ). The large sky coverage of the survey (almost 12,000 
deg 2 of sky) has resulted in phot o metri c measurements of approximately 357 



million objects (jAbazajian et al.l . 120091 ) . The completeness of SDSS catalogs 



for point sources is ~99.3% at the bright end and drops to 95% at magnitudes 
of 22.1, 22.4, 22.1, 21.2, and 20.3 in u, g, r, i, and z, respectively. Astrometric 
positions are accurate to better than 0. 1 arc-second per coordinate (rms) for 
sources with r < 20.5 ( jPier et al.l . 120031 ). and the morphologic al information 
from i mages allows re l iable star-galaxy separation to r ~ 21.5 (ILupton et al. 



2002; Scranton et al. 



20021 ). A compendium of the technical details about 
SDSS can be found on the SDSS web site ( jhttp: / /www.sdss.org[ ), which also 
provides the interface for public data access. 
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2.2 SDSS Photometric Measurements 



When discussing active (resolved) comets, the magnitudes used here are the 
"model magnitudes" measured by the SDSS. These are designed for galaxy 
photometry, and are determined by accepting the better of a de Vaucouleurs 



and a n exponential profile of arbitrary size and orientation (see lStoughton et al 



20021 ). W hile the surfac e brigh tness profile of a comet differs from that of 
a galaxy, ISolontoi et al.l (120101 ) demonstrate that these models still produce 
good fits to the cometary data. For the unresolved comets, the magnitudes are 
based on their point spread function (PSF) magnitude which is measured by 
fitting the point spread function model to the object. These two magnitudes 
serve as the basis for discriminating between resolved (galaxy-like) and unre- 
solved (star-like) sources in the SDSS; if the difference in magnitude between 
the r band PSF magnitude and the r band model magnitude is greater than 
0.145 then SPSS classifi es the object as resolved, or as a "galaxy" type object 



flStoughton et all 120021 ) . 



2.3 Finding Moving Objects in the SDSS Data 



Although mainly designed for observations of extragalactic sources, the SDSS 
is significantly contributing to studies of the solar system, notably in the 
success it has had with asteroid dete ctions, cataloged i n the SDSS Moving 



Object Catalog (hereafter SDSS MOC. Ilvezic et al.l . 120011 ). This public, value- 



added, catalog of SDSS asteroid observations contains, as of its fourth release, 
471,000 measurements of moving objects, 220,000 of which have been matched 
to 104,000 known asteroids from the ASTORB fikfU. The SDSS MOC data is 



Ivezic et al. 


i 


2001; 


Juric et al.. 


2002; 


Binzel et al.l. 


Assandri & Gil-Hutton, 


20081: 


Carvano et al.. 


2010) 



sec 



2007; Parker et al. 2008 



The SDSS camera ( IGunn et al.l . Il998f ) uses a drift-scan-like technique along 
great circles and detects objects in the order r, i, u, z, g, with detection in 
two successive bands separated in time by 72 seconds, with the diff erent bands 



regist ered to the same coordinate system using stationary stars (jPier et al. 



20031 ). Moving objects appear to have their photocenters spacially separated 
in different filters when color composite images are made. An angular velocity, 
calculated from the astrometrically-calibrated image centroids in each filter, 
is calculated for every photometric object. 

One of the features that allows the automatic creation of the SDSS MOC 
is that asteroids appear as point sources. A similar approach to the SDSS 



see ft p : / /ftp . lowell . edu/pub/ elgb / astor b . html 
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observations of resolved objects, in o rder to find comets, results in a sample 
well in excess of a million candidates (jSolontoi et al.l . 120101 ) , vastly larger than 
even the most optimistic estimate of observable comets. Even if every known 
comet in the sky were to be imaged, this sample would still be dominated by 
false positives by a factor of about 1000:1. Therefore, extended objects pose 
a much more difficult challenge for an automated reduction pipeline, further 
complicated by the fact that we are searching for moving extended objects, 
which to the SDSS photometric pipeline are interpreted as "moving galaxies." 



2.4 Finding Comets in the SDSS Data 



In order to successfully acquire comet observations from the SDSS database 
we employed two methods. The first involved making selection cuts based on 
SDSS measured photometry and data processing quality flags. The resulting 
candidate objects were then visually inspected (this method is discussed in 
depth by ISolontoi et al.ll2010f ). The advantage of this technique is that it is 
blind to the known comet sample. It does not matter whether the candidates 
are known comets, thought to be asteroids, or being observed for the first time 
(i.e. being discovered). The weakness of this approach is that it requires visual 
identification of each candidate object, and is subject to SDSS pipeline issues 
(e.g. deblending err ors, multiple epochs represented by a single image; see 
Solontoi et al.ll2010l for details). Since this technique was specifically designed 
for "cometary" (resolved) objects it will not select comets that are inactive, 
or have sufficiently low activity to appear as point sources at the resolution of 
the SDSS. 



A second technique t o select comets from the SDSS employs the method used 
by lJuric et al.l ( 120021 ) to id entify known aste roids in the SDSS MOC. We uti- 
lized the code developed by lJuric et al.l ( 120021 ) to propagate the orbits of known 
comets through the SDSS observational cadence. Much as in the case of the 
SDSS MOC, this code generates all possible RA and dec locations that a 
known comet could have been observed by the telescope over the course of the 
survey. This method enabled the selection of both resolved and unresolved 
comets. Visual inspection of candidates is still required as this list of positions 
sometimes has large astrometric uncertainties. Comets, unlike asteroids, have 
significant non-gravitational effects on their motion due to out-gassing during 
perihelion passage that alter the orbit of a comet over a given apparition, and 
certainly from passage to passage. For example, Comet Encke returns on av- 
erage about 2.5 hours sooner than predicted by gravitational com putations , 
a ph enomenon that was noticed as far back as the 19th century (IWhippld . 
1950h . and subsequent studies have demonstrated tha t the majority o f peri- 
odic comets experience non-gravitational acceleration (jSekaninal . If 9681 ). This 
fact, coupled with many comets not having well calibrated observations over 
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several orbits, leads to positional uncertainties in the predicted position. The 
comet can be too far from the predicted position to be uniquely matched by 
the algorithm, and in some cases it may not even be present in the SDSS ob- 
servation fi eld. A second weakness is that the underlying orbital code, based 



on OrbFit (jMilanil . Il999l ) can only be used to find comets with bound (ellip- 
tical) orbits, and thus excludes Long Period Comets. The final (and obvious) 
limitation is that it will only select known comets. 

Utilizing these two methods; 35 observations of 31 comets have been identified 
in the SDSS dataset, including two comets discovered by the survey (C/1999 
F2 Dalcanton and C/2000 QJ46 LINEAR). Two additional new comets were 
observed, but have not been successfully matched to any solar system body 
with a determined orbit, and thus will not be discussed in this paper; pho- 
tometry of those two objects is described by ISolontoi et al. fl201(j ). 



While 35 observations of 31 comets seems small in number compared the 
typical large number of objects assoc iated with the SDSS , they are in line 



with what is expected. As discussed in lSolontoi et al.l ( 120101 ). SDSS should be 



thought of as a single epoch survey in terms of comet results. While different 
comets are observed on different nights, the sky coverage for a single night is 
fairly small compared the entire sky, and the cadence of scans is not optimized 
to follow individual solar system objects. The usable area of the survey for 
comet work is about 20% of the whole sky. How many comets should we expect 
to find in such a survey? The Minor Planet Centej 2 1 lists approximately 200 
comets bright enough for robust detection (r < 20) on any given night. Taking 
these 200 potentially observable comets as a typical snapshot of the sky there 
should be approximately 40 comets in the survey (200 comets o n the sky * 20% 



of th e sky). Using SDSS DR5 to estimate our completeness (ISolontoi et al 



20101 ) we found it to be about 80%. Taking the sky coverage of SDSS and 
out completeness into account we expect to find approximately 32 comets, a 
prediction consistent with 35 observations of 31 comets reported here 

Table 1 lists the SDSS r band magnitudes, the difference between the measured 
PSF and model magnitudes (a point source would be equal to 0) and observing 
geometry (geocentric and heliocentric distances, and phase angle) for each of 
the comets discussed in this work. 



3 Analysis 



Here we discuss the methods for determining the photometric colors and sur- 
face brightness profiles for the resolved comets. From these measurements, 

2 http:/ /www. minorplanetcenter.org/iau/mpc. html 
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estimates can be made on the up per limit for nuclear r adius, and of the dust 
production via the quantity Afp (lA'Hearn et al.l . Il984j ). We then discuss the 
measurement of photometric colors and radii for the unresolved comets in 
section 3.2. 



3.1 Resolved Comets 



In this sample, 31 observations of 27 comets are classified by the SDSS as 
resolved sources, having a r band PSF— model magnitude value greater than 
0.145. For each of these comets the corrected frames (fits files) were obtained 
from the SDSS Data Archive Served" 3 ! (DAS) for all five filters. Nearby stars 
and galaxies were masked and radial surface brightness profiles were extracted 
independently in all five bands for each comet using routines in the IDP30 im- 
age manipulation package. These profiles consist of azimuthally averaged con- 
centric annuli about the optocenter of the comet in each photometric band. 
This allows four colors (u — g, g — r, r — i, and % — z) to be determined 
for all comets, even those whose SDSS assigned magnitudes are question- 
able due to photometric proc essing errors and non -optimal intensity profile 



for the flux measurement (see Solontoi et al.. 2010, sections 2.3, 3.4 and ref- 



erences therein). In the case of comets with well measured photometry, the 
total integrated magnitudes determined matched those reported by the SDSS 
photometric pipeline. 

As common with broad band filter studies of comets, we interpret the measured 
flux and resulting colors as due to the dust coma. Active comets are not 
free of gas however, and potentially strong emission lines may be present 
particularly in the wavelength ranges of the u and g bands. To verify these 
gas emission lines make a negligible contri bution to the integ rated flux (and 



therefore color) a sample comet spectru m (IHarris et all 120061) was convolved 



with the SDSS filter transmission curves (jDoi et all 120101 ) to produce u, g, and 
r— band magnitudes. Compared to a "smoothed" version of the spectrum (one 
with the emission lines removed) individual magnitudes changed by 3%, 0.1%, 
and < 0.1% for u, g, and r— band magnitudes respectively. These changes less 
than the standard deviation of the measured model magnitudes for the comets 
in each of these filters. 



3.1.1 Photometric Colors of Resolved Comets 



Extending on the work of ISolontoi et al.l ( 120101 ) we present the SDSS u — g 



g — r, r 



and i — z colors for the active comets observations in Table 2. 



3 http://das.sdss.org/www/html/ 

http://mips.as.arizona.edu/MIPS/IDP3 
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Figures 1 and 2 illustrate the uniformity in the colors of the active comets. 
Active comets occupy a very narrow distribution in color space, exhibiting the 
median colors of u — g : 1.57 ± 0.21, g — r : 0.57 ± 0.05, r — i : 0.22 ± 0.07, and 
i — z : 0.09 ±0.07. To aid in comparison to other Solar System bodies observed 
in the SDSS, we also calculate the SDSS asteroidal princi pal component colo r 



l a pc " defined as a pc 



0.89(# - r) + 0.45(r - i) - 0.57 (llvezic et all |200l|); 



the median a pc is 0.04 ± 0.06. 



This value lies between the modes for the a 



pc 



color distribution of main-belt asteroids (dominated by S and C types), and 
within uncertainty equal to the median value measured for the Jovian Trojans 
flSzabo et alll2007h . 



When these colors are transformed into BVRI photometric colors (llvezic et al. 



20071 ) they are shown be consisten t with previous color measurements of active 



comets (see 
of B-V: 0.74 



Solontoi et all |2010j, and references within) yeilding mean colors 
V— R: 0.44, and R— I: 0.58. These colors of active comets are 
slightly bluer (less than 1 — a) in B— V and V— R, and redder (less than 2 



ci- 



in R-I (-0.09 in B-V, -0. 06 in V-R, and +0 12 in R-I) than the colors 
found for cometary nuclei by lLamy &: Tothl (120091 ) . 



3.1.2 Surface Brightness Profiles 



For a si mple steady-state coma , the surface brightness should follow a p 1 
relation (IJewitt fc Meechl . Il987l ). with p being the projected linear distance 



from the nucleus. A plot of surface brightness (mag arcsec -2 ) against the 
log of angular s i ze, sh ould therefore show a logarithmic gradient m = — 1. 
Jewitt fc Meechl (119871 ) observed this relation and showed that when the effects 
of radiation pressure are taken into account, the slope can become steeper, 
(m = —1.5). Generally most measurements of the sur face brightness profiles 



of comets tend to fall in the range of — 2 < in < — 1 (IJewitt fc Meechl . 11987 



Lowry et all Il999l . for example). For comets that show profiles steeper than 



m 



-1.5 additional factors such as grain fading, may play a role. 



The slopes of the observed active SDSS comets' surface brightness profiles were 
fit using a weighted least-squares method independently in all five SDSS bands; 
we report these values for all resolved comets in Table 3. Figure 3 shows the 
surface brightness profiles and fits in all five bands for comet 67P/Churyumov- 
Gerasimenko, all of which are fit by m ~ —1. In addition, the colors are 
constant over the extent of the comet's profile. 

To express our confidence level in the measured surface brightness profiles, we 
have divided our comets into two "quality groups" based on the value of the 
difference between their PSF and model magnitudes. Comets that have small 
PSF— model magnitude values are less well resolved, and the accuracy of the 
fits to these surface brightness profiles are not as robust as for better resolved 
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comets. Figure 4 shows the trend of the mean slope (the mean of the g, r and 
i slopes) as a function of PSF— model magnitude, and demonstrates a break 
near PSF-model = 1.5. At this point, the less resolved comets begin a trend 
toward steeper slopes that extends all the way down to point-sources (PSF- 
model ~ 0). These point sources (star symbols in Figure 4) are stars from 
the comet fields of similar magnitude and were fit by the same slope fitting 
routine used for the surface brightness profiles of the comets. Since it is not 
clear which of the comets are less resolved due to low activity versus those due 
to the resolution limit of SDSS, we will use a value of PSF— model magnitude 
difference of 1.5 to define two "quality groups" of surface brightness profile 
slopes. The well resolved comets (PSF-model > 1.5) are placed in QG1 and 
the less resolved comets (PSF— model < 1.5) in QG2 



For several comets, the u and z bands do not have high enough surface flux 
density to allow for robust fits to be made, and in those cases no fit is given in 
Table 3. The g,r, and i bands are generally self consistent for a given comet. 
This trend of similar slopes from the g, r, and i band fits is illustrated in Figure 
5. The solid lines trace a 1:1 relationship. The fits to the g, r, and i band slopes 
trace each other well, while the fits to the u and z bands show more scatter. 
The lower quality of fits in the u and z bands is due largely to the surface flux 
density not being large enough compared to the sky to allow for rigorous fits. 
This results from a combination of the comets themselves being fainter in the 
u band, and the lower q uantum efficiency of the SPS S camera for the u and 



qi 

z filters flFukueita et all . Il996t IStouehton et all . l2002h . 



These five-band surface brightness profiles reveal that not only are the pro- 
file slopes consistent across a wide range of wavelengths, as seen in Figure 
5 and Table 3, bu t changes in the profile slope are also consistent from one 
band to the next. IJewitt fc Meechl (119871 ) predicted a transition to a steeper 
profile slope at a distance from the nucleus where solar radiation pressure 
becomes dominant. Several bright comets observed by the SDSS demonstrate 
this behavior, as illustrated by comet 67P/Churyumov-Gerasimenko in Figure 
3. The slope transition is seen at the same distance (approximately 6 x 10 5 km 
in the case of 67P) for a given comet across all filters where the flux density 
is significantly above the background. 



3.1.3 Upper Limits to Cometary Radii 



Upper-limit estimates of the radius of the cometary nuclei may be determined 
by using the PSF fitted magnitude for the resolved comets. We use the mea- 



sured PSF magnitude from the SDSS photometric pipeline (IStoughton et al. 



20021 ). which involves sync-shifting the image so that it is exactly centered on 



a pixel, and then fitting a Gaussian model of the PSF to it, with additional 
corrections applied in order to take into account the full variation of the PSF 
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across the field. It would not be useful to take flux from a large aperture for 
these determinations, as the light from the coma is vastly more luminous than 
that from the nucleus. While the detected flux in the PSF-magnitude does 
indeed come from light interacting with the coma, the PSF-magnitude repre- 
sents the most light that could be theoretically coming from a point source in 
the region of the nucleus. In a sense this is the same assumption that is made 
during the calculation of radii from photometry of unresolved comets. An un- 
resolved comet does not show a coma at the level of photometric resolution, 
and so the PSF photometry of the object is taken to be the light reflecting 
from the nucleus, even though that light may in fact be from an unresolved 
coma, rather than the nucleus itself. 

Keeping in mind these assumptions, one can take the PSF-magnitude mea- 
surement in the r band as a nuclear upper limit, and calculate the absolute 
magnitude ?"(i,o) : the r band magnitude the comet would have with a geocen- 
tric and heliocentric distance of 1 AU and viewed at a phase angle of zero 
degrees (this is a physically impossible geometry). So for an observed r band 
magnitude 

r (1 , O )=r-5bg[flA]-0(a) (1) 



where R and A are the heliocentric and geocentric distances in AU, a is the 
observed phase angle, and 0(a) is the phase function. This absolute magnitude 
is constructed in the same way using the SDSS r band magnitude as is the 
standard IAU asteroid H magnitude with the Johnson V band magnitudes 
(for reference r^o) ~ H — 0.2). This absolute magnitude can be w ritten in 



terms of the comet's diameter, D and albedo, A (llvezic et all 120011 ). 



r (i,o) 



17.9 



2.51og(— ) 

8V o.r 



51og( 



1km' 



(2) 



By assuming a standard albedo of 0.04 and adopting a linear phase function, 
0(a) = (3a, with (3 = 0.035 mag deg -1 , these two equations can be used to 
transform the PSF r-band magnitudes into upper limits for the sizes of the 
comet nuclei. The upper radii limits derived in this manner (Table 4) are 
comparable to JFC radii in the lite rature, in parti c ular t he list of Jupiter 
Family Comet size estimates made by iTancredi et al.l ( 120061 ) , who also assume 
an albedo of 0.04. The t rend is that the SPSS upper limits tend to be larger 
than those published by ITancredi et al.l (120061) by about 75%, particularly for 
the comets for which ITancredi et al.l ( 120061 ) express high levels of confidence 
in their estimates. 



3.1.4 Dust Production Rates 

The dust production rate of comets gives insight into the current status, and 
if measured over time, the evolution of the comet's activity. Dust production 
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rates for the observed comets may b e characterized from their photometry 
through the quantity Afp defined by lA'Hearn et al.l (119841 ) . Where A is the 
dust grain albedo, / the filling factor within the chosen aperture, and p the 
linear radius corresponding to the aperture. Afp can be computed directly 
from observable quantities 



Afp(cm) 



(2AR) 
P 



2 p 

" comet 



T 



(3) 



Here A (cm) and R(AXJ) are the geocentric and heliocentric distances, p(cm) 
is the linear radius of the photometric aperture, and F comet /F Q is the ratio of 
the observed comet flux to that of the Sun. Empirical comparisons indicate a 
linear relationship between the Afp value and the dust production rate, with 
Afp (measured in units of 1000 cm) being roughly equal to the modeled dust 



production rate in metric tons per second ( lA'Hearn et al.l . Il995l ) 



For an idealized steady-state coma, the surface brightness profile follows a p -1 
profile. Under that coma model, Afp is independent of aperture, but since real 
comae deviate from m = — 1 it is necessary to define the radius at which Afp is 
to be measured. Rat her than a s sumin g an arbitrary radius we take a cue from 
galaxy photometry (jPetrosianl . Il976l ) and define the radius, p at which Afp 
is calculated as the radius at which the local surface brightness value is 10% 
of the mean enclosed local surface brightness value, m = 10%. This 

° ' mean[b B(<p)\ 

criterion results in Afp being measured at approximately 10 3 — 10 5 km from 
the center of the comet, and is consistent with radii chosen in the literature 



for comets over a wide range of o bservational distances (cf. lLowry et al.l 12003 
and iMazzotta Epifani et al.l 120071 ) . 



While under idealized conditions Afp is independent of aperture (for a 
profile), and wavelength (for a grey dust approximation), it is dependent on 
the observational p hase angle. B y adopting a simple phase function for single 
dust particles from iDivind (Il98ll ). the measured Afp values may be corrected 
to A(0)/p, the expected value if the comet had been at zero phase angle dur- 
ing the observation. As most reported Afp data in literature are not phase 
corrected, we will generally refer to our uncorrected Afp values, unless specif- 
ically referencing the phase-adjusted A(0)/p. Table 4 lists the Afp and A(0)/p 
measurements for the active comets. 



3.2 Unresolved Comets 



Four comets in our sample are unresolved by SDSS. Two of these, 19P/Borrelly 

and 113P/Spitaler, are dynamically classified as Jupiter Family Comet s. 174P/Ech eclus 

is the Centaur 60558 Echeclus that has been observed to have activity ( Choi et all . 
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20061 ). and 176P /LINEAR is a main belt comet, previo usly designated as main 



belt asteroid 118401 LINEAR (IHsieh &: Jewittl . 120061 ) . All four of these unre- 
solved comets are in the SDSS MOC; 124P and 126P are correctly identified 
and matched to their asteroid designation, while 19P and 113P, which are in 
comet but not asteroid databases, are unmatched sources in the SDSS MOC. 
There are likely to be more unresolved comets to be found in the SDSS MOC, 
but identifying and retrieving these objects is difficult due to non-gravitational 
forces causing them to deviate from projected positions on the sky. 



The u — g, g — r, r — i, and i — z colors are presented in Table 2. Figure 2 shows 
that the unresolved JFCs, 19P/Borrelly and 113P/Spitaler, fall into the color 
space spanned by the active comets. 

If we assume that these observations of unresolved comets are images of the 
inactive cometary nucleus, we can make an estimate of the comet's size as- 
suming an albedo of 0.04. These are included in Table 4, and like the upper 
limits estimated for the resolved comets, are in satisfactory agreement with 
previously determined values. 



4 Discussion 



4-1 Photometric Colors 



In light of the significant differences of observational geometry, physical pa- 
rameters, and orbital type, all active comets have the same SDSS colors, which 
span a range of ~ 6000 A across the five filters. The photometric (u — g, g — r, 
r — i, and i — z) colors of the active comets (Table 2) reveal that they are 
remarkably similar to each other, particularly in g — r, r — i, and i — z, where 
the colors have a scatter of < 0.07 magnitude. Furthermore, this uniformity of 
color does not show systematic variation with observational geometry (phase 
angle, geocentric or heliocentric distance), measured quantities (r band magni- 
tude, PSF— model magnitude, slope of the surface brightness profile), derived 
quantities (nuclear radius, A/p), or orbital properties (orbital inclination, ec- 
centricity and perihelion distance). Figure 6 illustrates the dependence of all 
four colors on each of these parameters. The colors are consistent with no slope 
within la against all parameters as shown in Figure 7. The result of the colors 
of active comets ha ying no correlati o n wit h heliocentric distance is consistent 
with the results of IJewitt fc Meechl (119881 ) . Considering that the dust in the 
coma is thought to be primordial solar system grains freed by the volatiliza- 
tion of ices, this result may indicate that such a color uniformity could be a 
common property of the primordial dust grains of the outer solar system. 
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The colors of the comets discussed h ere are strikingly similar to that of the 
Trojans of Jupiter ( jSzabo et all 120071 ) as shown in Figure 2. Trojans al s o show 
similar colors to active comets in 5.2-38/zm observations lEmery et al.l ( 120061 ). 
which may indicate a link between the surface chemistry and composition of 
comets and Trojans. 



To further characterize this color similarity we examine the gradient of the 
normalized reflectivity, S , expressed as % per 1000A. The colors and dis- 
tribution of S for the active comets are compared to those of a sample of 
363 Trojans, 12 Centaurs, and 23 Trans Neptunian Objects (TNOs) from the 
SDSS. The Trojan and Centaur samples were selected by their orbital elements 
(and positive cross-identific ation) from the SD SS MOC, and the TNOs were 
found in SDSS Stripe 82 by lBecker et all (I201CH ) . Figure 8 shows the range in 
color-color space occupied by these populations, with the Trojans and active 
comets being in close agreement . The SDSS colors found in this work agree 
with those of ISzabo et al.l ( 120071 ). and the S for t he act ive comets resembles 
the S distribution of JFC nuclei found by I Jewittl ( 120021 ). Note that Centaurs 
are extremely faint objects with correspondingly higher photometric errors. 



Histograms of these three populations are seen in Figure 9. The faintness of 
the Centaurs and TNOs makes the full five-band gradient less reliable. The 
colors and reflectivity gradients for these populations are summarized in Ta- 
ble 5. By visual inspection, these histograms confirm in S what was seen in 
their colors: the active comets and Trojans are remarkably similar. To char- 
acterize this similarity a Kolmogorov-Smirnov test (KS-test) was performed 
on the S distributions of the populations in question. Between the Trojans 
and the Comets, the KS-test probability value was 0.16 with respect to the S 
distribution across all five bands, and rises to 0.81 for the gri band, suggest- 
ing that the comets and the Trojans could have been drawn from the same 
population. To reiterate, even though Trojans and Comets are two different 
types of object, with the light we see reflecting off of different sources (solid 
body vs. dusty coma), their colors are statistically indistinguishable. 



4-2 Nuclear Radii 



Estimates of absolute magnitude and nuclear size were made using equations 
1 and 2. Despite the measurements of nuclear radii for the active comets being 
the result of accepting an upper limit from the SDSS PSF magnitude, the radii 
presented here agree with the measurements in literature as shown in Table 
4. 

Even for targeted comet surveys, it is difficult to be certain that all of the flux 
seen is really light that has interacted with the surface of the comet; in some 
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cases there may be low-level activity below the SDSS resolution limit. Cer- 
tainly there are methods of obtaining highly accurate radius measurements, 
the best obviously b eing in si t u spa cecraft encounte rs with comets. Only fiv e 
comets, lP/ Hallev (iReinhardl. llQp) , 19P/Borrell v dSoderblom et all l2002t ). 
81P/Wild-2 flBrownlee et al.l. l2004h. 9P/ Tempel-1 flA'Hearn et"aD . l2005h . and 
103P /Hartley- 2 flA'Hearn & DIXI Teaml . l201lh have to date had observations 
of the nucleus by remote spacecraft. 



Lamy fc Tothl (119951 ) outline how to use the high spatial resolution afforded by 
the Hubble Space Telescope to disentangle the nucleus signal from the inner 
coma. Highly accurate results may also be produced if one has access to si- 
multaneous observa tions of the comet at both visible and infrared wavelengths 
(jLamy et all 120021 ). Most cometary observations, however, must make use of 
similar methods as employed here. For resolved comets the coma must be 
algorithmically removed, often either through the modeling of the coma and 
subtraction, or through the use of small aperture photometry or PSF model 
fitting. Even unresolved comets must be assumed to be inactive for observa- 
tion of the radius. Complicating the matter further is the non-uniformity of 
cometary albedos: while comets are dark overall, the albe do can vary comet-to - 
comet, and even across the surface of an individual one ( jBuratti et al.l . 120041 ) . 
Being small collisionally-shaped bodies, they can have axial, and even tri-axial 
shapes. This feature presents an even larger problem as very few comets have 
had their nuclear light curve accurately measured, (only a few dozen comets 
have robustly measured nucleus observations) and many of those show dra- 
matic asymme t ries ( e.g. comet 19P/Borrelly with a primary axis ratio of 2.5:1 



flBuratti et all 120041 ) I 



Although the SDSS resolution is not sufficient for direct measures of the radii 
they do have the advantage that all observations are made and reduced by a 
well calibrated and highly characterized process. Thus, the large number of 
comets observed by the SDSS allows a statistical study of the JFC popula- 
tion, with confidence that the data set is self consistent. To better compare 
these results with previous cometary studies, we have used the estimated radii 
upper limits of the comets to calculate (assuming an albedo of 0.04) the more 
commonly used absolute H magnitude. Figure 10 shows the Cumulative Lu- 
minosity Function (CLF - the number of JFCs more luminous than a given 
H), and the Cumulative Size Distribution (CSD) for the JFCs. These distri- 
butions may be fit with power laws, with the CLF best fit by an exponent of 
0.49 ± 0.05, and the CSD by -2.4 ± 0.2 for H < 18. These give the functional 
form of the distributions as 



CLF: A^(< H) oc 10 a49H 
CSD: N{< r) oc r' 2A 



(4) 
(5) 
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for a given absolute magnitude H and radius r. These values are in line with 
the results of previous studies, which are summarized in Table 6. 



We find that these distributions may also be fit with a broken power law, with 
breaks at H = 14.5 in the CLF, and at radius = 4km in the CSD (Figure 10). 
For the CLF the gradients are 0.73 ± 0.08 for H < 14.5 and 0.19 ± 0.03 for 
14.5 < H < 18 yielding -3.1 ±0.4 and -1.0 ±0.1 for the CSD. This is similar 



to the CLF measured by lLowry et al.l ( 120031 ) when considering their entire 
sample (both resolved and unresolved) of comets (0.53±0.04 and 0.22±0.02). 
We find that our shallow power law for smaller comets are is agreement, but 
find that our sample supports a steeper fit for the larger comets. Again these 
slopes are derived from both radii found for unresolved comets, and upper 
limits placed on the nuclear radii of resolved comets. While the parameters of 
this broken power law fit are not strongly constrained they do suggest that 
that faint comet population could be much shallower than indicated by a single 
power law fit. 



Table 7 compares our results for the Jupiter Family Comets with those of 
other small body populations. As is the case for many of these p opulations, 
the ob served CLF power law exponent is similar to 0.5, predicted by iDohnanyi 
( 119691 ) for a model based on an equilibrium cascade of self-similar collisions. 
The large number of Main Belt Asteroids observed by the SDSS allows their 
CLF to be analyzed family-by-family, and an alytically fit to a fu nction that 
accounts for a change in the power law slope (IParker et al.l . 120081 ) . With the 
multitude of observations on other small body populations expe cted from the 
next generation o f sky s urveys (e.g. Pan-STARRS and LSST, Kaiser et al 



20021 ; llvezic et al.l . 120081 ). it is expected that similar detailed analysis could 



be performed on the less well-sampled populations including Jupiter Family 
Comets. 



4.3 Afp 



The Afp values obtained for the active comets observed by the SDSS cover a 
range of three orders of magnitude 0.88 cm < Afp < 817.68 cm. Table 8 com- 
pares the Afp values determined in this paper with those from other studies. 
They are consistent with other measurements; in particular 67P/Churyumov- 
Gerasimenko (the ESA Rosetta Mission target) has a well sampled Afp over 
the comet's orbit. Our A fp value for this comet is consistent with that found 
at the same phase angle by Agarwal et al. ( 2010l ). 



Using Afp as proxy for dust production, these values show a trend of the 
non-JFC comets having significantly higher dust production than do JFC 
comets. This supports the suggestion that Long Period Comets being recently 
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introduced to the inner solar system, are not as surf ace- volatile depleted as a 
Jupiter Family Comet that would have presumably made multiple perihelion 
passages (cf Figure 11). Other factors such as a difference in the size distribu- 
tion or albedo of the dust between these two populations could also produce 
a change in Afp. 



5 Summary 



We have analyzed data for 31 comets observed in the SDSS, which provides ac- 
curately calibrated and measured multi-wavelength photometry in five bands. 
These comets span a wide range of heliocentric distances, observational, and 
orbital parameters and all have been observed with the same instrument and 
processed with the same software. This photometry has been used to make 
measurements of the colors, sizes, surface brightness profiles and rates of dust 
production (in terms of the Afp formalism). Our main results are as follows. 

(1) Despite the variety of cometary parameters, the distributions of photo- 
metric colors are extremely narrow (0.57 ±0.05 in g — r for example), and 
statistically indistinguishable from those of Jupiter Trojans. The comets 
exhibit no correlation between color and physical, dynamical, or observa- 
tional parameters, as seen in Figure 7. Additionally the surface brightness 
profile for each comet is found to be invariant with wavelength in the op- 
tical. 

(2) The uniform, red color of the comets indicates that the light from the 
dust in the coma is complex, no showing any systematic variation with 
observational geometry, dynamics or physical properties of the comets, 
and can not be explained through simple scattering (e.g. a simple 1/A n 
optical depth dependance with n ~ 4). The visible light observed here is 
light scattered by the larger particles in the coma, either conglomerates of 
small grains, or larger, macroscopic rocks. The similarity in reflectance to 
solid surfaces (both cometary nuclei, and Jupiter Trojans) may indicate 
that the material responsible for the dark, reddish appearance of these 
solid objects is also present in the coma and dictates its optical properties. 

(3) We find that the cumulative luminosity function for the Jupiter Family 
Comets can be fit by a power law of the form iV(< H) oc io(°- 49±a05 )^, 
with evidence of a broken power law with a exponent of 0.73 ± 0.08, 
transitioning to 0.19 ± 0.03 at H ~ 14.5, suggesting a shallower gradient 
for the faint population. These results are consistent with distributions 
of small bodies, both of JFCs, and other Solar System populations. 

The analysis presented here is relevant to upcoming; large- sky surveys such 



as th e Dark Energy Survey (IFlaugher fc Dark Energy Survey Collaboration 



20071 ). Pan-STARRS (IKaiser et al.l . 120021 ) and the Large Synoptic Survey Tele 
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scope (jlvezic et all 12008k LSST). The results of this work represent a "snap- 
shot" survey of come t s, obs erving a sample of comets at a single point in their 
orbit. Solontoi et al. (bold ) estimate that 10 3 to 10 4 comets will be observed 
in a survey such as LSST, with a limiting magnitude of r ~ 24.5. Further, 
over the course of 10 years LSST will observe each come t repeatedly. Base d 



on the results of a simulated LSST observational cadence (jlvezic et all 120081 ) 



a JFC such as 70P/Kojima would be observed 300-400 times, spanning its 
whole heliocentric range. Such observations mean that not only may the anal- 
ysis presented here be done for a vast number of comets, down to even fainter 
limiting magnitudes, but the time evolution of each of these parameters may 
be studied over the comet's orbit. 
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0.01 


1.90 


3.559 


4 


338 


8.9 


67P/Churyumov-Gcrasimcnko 


14.28 




0.01 


3.15 


1.590 


1 


836 


32.5 


69P/Taylor 


15.59 


± 


0.01 


3.19 


1.127 


1 


950 


21.0 


70P/Kojima 


16.68 




0.01 


1.97 


1.629 


2 


594 


7.6 


tll3P/Spitaler 


21.05 




0.05 




2.332 


3 


320 


2.3 


129P/Shoemaker-Lovy 


18.07 




0.02 


2.03 


2.495 


3 


051 


17.0 


146P/Shocmakcr-Linear 


18.78 




0.02 


1.27 


1.301 


2 


023 


23.9 


158P/Kowal-LINEAR (1) 


18. 78 


i 


0.02 


1.15 


3.619 


4 


596 


3.0 


158P/Kowal-LINEAR (2) 


18.60 


i 


0.02 


1.10 


3.650 


4 


597 


4.2 


158P/K 


owal-LINEAR (3) 


19.36 


± 


0.02 


1.01 


3.988 


4 


798 


7.3 


fl74P (60558 Echoclus) 


21.20 


± 


0.08 




14.350 


IE 


.167 


2.27 


fl76P (118401 LINEAR) 


20.27 


± 


0.03 




2.374 


3 


256 


9.2 


P/2002 E.I57 (LINEAR) 


18.99 


± 


0.01 


0.37 


1.752 


2 


720 


5.9 


2003 WY25 


18.33 


± 


0.01 


0.66 


0.265 


1 


217 


25.8 


C/1999 


F2 (Dalcanton) 


15.81 


± 


0.01 


2.29 


4.344 


4 


996 


9.3 


C/2000 


K2 (LINEAR) 


16.82 


± 


0.01 


1.74 


3.950 


4 


776 


7.2 


C/2000 


QJ46 (LINEAR) (1) 


17.44 


± 


0.01 


1.64 


1.162 


2 


168 


2.3 


C/2000 


QJ46 (LINEAR) (2) 


19.39 


± 


0.04 


1.47 


3.318 


3 


657 


15.3 


C/2000 


SV74 (LINEAR) 


14.68 


± 


0.00 


3.55 


4.015 


4 


380 


12.6 


C/2000 


Y2 (Skiff) 


16.71 


± 


0.01 


1.19 


1.824 


2 


785 


5.8 


C/2001 


RX14 (LINEAR) 


12.62 


± 


0.01 


2.76 


1.619 


2 


103 


26.9 


C/2002 


07 (LINEAR) (1) 


19.37 


± 


0.03 


1.23 


5.645 


5 


886 


9.8 


C/2002 


07 (LINEAR) (2) 


15.24 


± 


0.00 


2.94 


2.504 


3 


041 


17.4 


P/1999 


VI (Catalina) 


17.36 


± 


0.01 


1.87 


2.164 


3 


099 


7.3 


P/2002 


T5 (LINEAR) 


18.33 


± 


0.02 


1.97 


4.351 


5 


002 


9.1 


P/2004 


Al (LONEOS) 


18.55 


± 


0.01 


0.99 


5.226 


5 


490 


10.4 


P/2006 


U5 (Christensen) 


17.63 


± 


0.01 


2.04 


1.538 


2 


358 


16.6 



Table 1 

The observed comets. The r band magnitudes reported are SDSS model magnitudes 
except in the case of unresolved comets (denoted by a f) where the SDSS PSF 
magnitudes are reported. A, R, and phase list the Geocentric and Heliocentric 
distance (AU) and phase angle of the comet at the time of observation. 
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Comet 


r 


u 




9 


g — r 


r — i 


i — z 


ACTIVE COMETS 


















30P Rcinmuth 


14.85 ± 


0.01 


1.57 ± 0.02 


0.60 ± 0.01 


0.32 ± 0.01 


0.19 ± 0.01 


46P/Wirtancn 


18.84 ± 


0.02 








0.54 ± 0.09 


0.25 ± 0.05 


0.06 ± 0.05 


47P /Ashbrook- Jackson 


20.57 ± 


0.05 


1.09 ± 0.33 


0.59 ± 0.07 


0.38 ± 0.07 


0.02 ± 0.17 


50P/Arcnd 


18.13 ± 


0.03 


1.72 ± 0.32 


0.57 ± 0.04 


0.08 ± 0.04 


0.27 ± 0.10 


62P/Tsuchinshan 


15.04 ± 


0.01 


1.27 ± 0.04 


0.45 ± 0.01 


0.18 ± 0.01 


0.07 ± 0.02 


64P/Swift-Gehrels 


19.23 ± 


0.04 


1.8 ± 0.45 


0.60 ± 0.06 


0.22 ± 0.06 


0.00 ± 0.21 


65P/Gunn (1) 


17.13 ± 


0.01 


1.68 




0.07 


0.58 ± 0.01 


0.21 ± 0.01 


0.08 ± 0.03 


65P/Gunn (2) 


17.15 ± 


0.01 


1.57 


± 


0.06 


0.55 ± 0.01 


0.21 ± 0.01 


0.10 ± 0.03 


67P/Churyumov-Gcrasimcnko 


14.28 ± 


0.01 


1.83 


± 


0.03 


0.70 ± 0.01 


0.25 ± 0.01 


0.06 ± 0.01 


69P/Taylor 


15.59 ± 


0.01 


1.51 


± 


0.04 


0.60 ± 0.01 


0.25 ± 0.01 


0.06 ± 0.02 


70P/Kojima 


16.68 ± 


0.01 


1.47 


± 


0.07 


0.64 ± 0.01 


0.23 ± 0.01 


0.21 ± 0.03 


129P/Shocmakcr-Levy 


18.07 ± 


0.02 


1.18 


± 


0.12 


0.55 ± 0.03 


0.17 ± 0.03 


0.09 ± 0.07 


146P/Shocmakcr-Lincar 


18.78 ± 


0.02 


1.46 


± 


0.15 


0.64 ± 0.03 


0.25 ± 0.03 


0.12 ± 0.06 


158P/Kowal-LINEAR (1) 


18.78 ± 


0.02 


1.59 


± 


0.19 


0.56 ± 0.03 


0.22 ± 0.03 


0.01 ± 0.09 


158P/Kowal-LINEAR (2) 


18.60 ± 


0.02 


1.48 


± 


0.14 


0.59 ± 0.03 


0.19 ± 0.02 


0.20 ± 0.05 


15SP/Kowal-LINEAR (3) 


19.36 ± 


0.02 


2.09 


± 


0.39 


0.58 ± 0.03 


0.20 ± 0.04 


0.19 ± 0.08 


P/2002 EJ57 (LINEAR) 


18.99 ± 


0.01 


1.57 


± 


0.11 


0.54 ± 0.02 


0.20 ± 0.02 


0.12 ± 0.05 


2003 WY25 


18.33 ± 


0.01 


1.47 


± 


0.08 


0.64 ± 0.01 


0.27 ± 0.01 


0.12 ± 0.04 


C/1999 


F2 (Dalcanton) 


15.81 ± 


0.01 


1.69 


± 


0.05 


0.52 ± 0.07 


0.25 ± 0.05 


0.07 ± 0.05 


C/2000 


K2 (LINEAR) 


16.82 ± 


0.01 


1.55 


± 


0.06 


0.55 ± 0.01 


0.22 ± 0.01 


0.10 ± 0.02 


C/2000 


QJ46 (LINEAR) (1) 


17.44 ± 


0.01 


1.55 




0.11 


0.57 ± 0.02 


0.23 ± 0.01 


0.12 ± 0.04 


C/2000 


QJ46 (LINEAR) (2) 


19.39 ± 


0.04 


1.04 


± 


0.27 


0.56 ± 0.06 


0.41 ± 0.06 


0.09 ± 0.13 


C/2000 


SV74 (LINEAR) 


14.68 ± 


0.00 


1.66 


± 


0.02 


0.52 ± 0.01 


0.21 ± 0.01 


0.04 ± 0.01 


C/2000 


Y2 (Skiff) 


16.71 ± 


0.01 


1.58 


± 


0.04 


0.55 ± 0.01 


0.18 ± 0.01 


0.09 ± 0.01 


C/2001 


RX14 (LINEAR) 


12.62 ± 


0.01 


1.61 


± 


0.03 


0.57 ± 0.02 


0.23 ± 0.02 


0.06 ± 0.01 


C/2002 


07 (LINEAR) (1) 


19.37 ± 


0.03 


1.35 


± 


0.21 


0.58 ± 0.04 


0.14 ± 0.04 


-0.12 ± 0.25 


C/2002 


07 (LINEAR) (2) 


15.24 ± 


0.01 


1.56 


± 


0.07 


0.56 ± 0.02 


0.22 ± 0.02 


0.07 ± 0.03 


P/1999 


VI (Catalina) 


17.36 ± 


0.01 


1.73 


± 


0.09 


0.58 ± 0.01 


0.23 ± 0.01 


0.09 ± 0.03 


P/2002 


T5 (LINEAR) 


18.33 ± 


0.02 


1.52 


± 


0.21 


0.66 ± 0.03 


0.23 ± 0.03 


0.09 ± 0.10 


P/2004 


Al (LONEOS) 


18.55 ± 


0.01 


1.67 


± 


0.20 


0.61 ± 0.02 


0.31 ± 0.02 


0.09 ± 0.05 


P/2006 


U5 (Christcnsen) 


17.63 ± 


0.01 


1.46 


± 


0.09 


0.51 ± 0.01 


0.08 ± 0.01 


0.22 ± 0.04 


Median Color 






1.57 ± 0.21 


0.57 ± 0.05 


0.22 ± 0.07 


0.09 ± 0.07 


UNRESOLVED COMETS 


















19P/Borrelly 


21.35 ± 0.05 


1.11 ± 0.50 


0.77 ± 0.09 


0.12 ± 0.09 


0.27 ± 0.25 


113P/Spitaler 


21.05 ± 0.05 


1.06 ± 0.44 


0.60 ± 0.08 


0.33 ± 0.07 


0.00 ± 0.61 


174P (60558 Echcclus) 


21.20 ± 0.08 


1.55 


± 


1.02 


0.73 ± 0.13 


0.18 ± 0.08 


-0.71 ± 0.53 


176P (118401 LINEAR) 


20.27 ± 0.03 


1.67 


± 0.40 


0.54 ± 0.04 


0.06 ± 0.05 


0.04 ± 0.10 



Table 2 

The SDSS colors of the comets. The listed r band magnitudes are SDSS model 
magnitudes for active comets, and SDSS PSF magnitudes for the unresolved comets. 
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Comet 


Psf-model 


u 


Surface 

a 


Brightness Profi 
r 


le Slope 

i 


2 


30P Re 


nmuth 


2.20 


-1.57 


± 


0.16 


-1 


57 


± 


0.01 


-1.58 


± 


0.01 


-1.54 


± 


0.01 


-1.55 


± 


0.01 


46P/Wirtancn 


1.63 




- 




-1 


18 


± 


0.06 


-1.22 


± 


0.06 


-1.34 


± 


0.06 


-1.54 


± 


0.40 


47P/ Ashbrook- Jackson 


0.24 




- 




-1 


87 


± 


0.26 


-2.85 


± 


0.28 


-1.74 


± 


0.06 


-2.43 


± 


0.31 


50P/Arcnd 


2.35 




- 




-1 


31 


± 


0.07 


-1.38 


± 


0.11 


-1.50 


± 


0.15 


-1.32 


± 


0.29 


62P/Tsuchinshan 


3.91 


-0.91 


± 


0.07 


-0 


86 


± 


0.01 


-0.97 


± 


0.01 


-0.98 


± 


0.01 


-0.90 


± 


0.03 


64P/Swift-Gehrcls 


1.75 




- 




-1 


29 


± 


0.18 


-1.22 


± 


0.11 


-1.28 


± 


0.16 




- 




65P/Gunn (1) 


1.75 


-1.58 


± 


0.10 


-1 


88 


± 


0.05 


-1.94 


± 


0.04 


-1.94 


± 


0.04 


-1.66 


± 


0.17 


65P/Gunn (2) 


1.90 


-2.27 


± 


0.56 


-1 


87 


± 


0.03 


-1.91 


± 


0.04 


-1.82 


± 


0.04 


-1.73 


± 


0.07 


67P/Churyumov-Gcrasimcnko 


3.15 


-0.99 


± 


0.04 


-1 


05 


± 


0.01 


-1.04 


± 


0.01 


-1.02 


± 


0.01 


-1.05 


± 


0.02 


69P/Taylor 


3.19 


-1.38 


± 


0.24 


-1 


24 


± 


0.01 


-1.35 


± 


0.01 


-1.33 


± 


0.01 


-1.42 


± 


0.04 


70P/Kojima 


1.97 


-0.97 


± 


0.41 


-1 


32 


± 


0.03 


-1.41 


± 


0.02 


-1.43 


± 


0.02 


-1.39 


± 


0.06 


129P/Shoemaker-Lcvy 


2.03 








-1 


85 


± 


0.10 


-1.67 


± 


0.04 


-1.53 


± 


0.07 


-1.13 


± 


0.20 


146P/Shoemaker-LINEAR 


1.27 








-1 


62 


± 


0.12 


-1.74 


± 


0.11 


-1.74 


± 


0.10 


-1.75 


± 


0.42 


158P/Kowal-LINEAR (1) 


1.15 








-1 


99 


± 


0.11 


-2.20 


± 


0.10 


-2.03 


± 


0.15 


-2.10 


± 


0.19 


158P/Kowal-LINEAR (2) 


1.10 


-2.28 


± 


0.59 


-2 


30 


± 


0.11 


-2.30 


± 


0.08 


-2.08 


± 


0.10 


-2.79 


± 


0.49 


15SP/Kowal-LINEAR (3) 


1.01 








-2 


12 


± 


0.17 


-1.85 


± 


0.07 


-1.83 


± 


0.11 


-1.77 


± 


0.33 


P/2002 EJ57 (LINEAR) 


0.37 


-1.87 


± 


0.36 


-2 


84 


± 


0.18 


-2.96 


± 


0.14 


-3.01 


± 


0.12 


-3.14 


± 


0.27 


2003 WY25 


0.66 


-1.24 


± 


0.39 


-1 


88 


± 


0.07 


-1.76 


± 


0.06 


-2.04 


± 


0.09 


-1.90 


± 


0.17 


C/1999 


F2 (Dalcanton) 


2.29 


-1.12 


± 


0.10 


-1 


08 


± 


0.01 


-1.07 


± 


0.01 


-1.03 


± 


0.01 


-1.07 


± 


0.03 


C/2000 


K2 (LINEAR) 


1.74 


-2.03 


± 


0.41 


-1 


94 


± 


0.03 


-1.92 


± 


0.03 


-1.89 


± 


0.02 


-1.71 


± 


0.05 


C/2000 


QJ46 (LINEAR) (1) 


1.64 


-2.05 


± 


0.35 


-1 


69 


± 


0.04 


-1.73 


± 


0.02 


-1.89 


± 


0.04 


-1.95 


± 


0.13 


C/2000 


QJ46 (LINEAR) (2) 


1.47 








-1 


96 


± 


0.32 


-1.49 


± 


0.21 


-1.96 


± 


0.44 


-1.96 


± 


0.60 


C/2000 


SV74 (LINEAR) 


3.55 


-1.19 


± 


0.17 


-1 


24 


± 


0.01 


-1.23 


± 


0.01 


-1.22 


± 


0.01 


-1.20 


± 


0.02 


C/2000 


Y2 (Skiff) 


1.19 


-2.08 


± 


0.22 


-2 


11 


± 


0.02 


-2.17 


± 


0.03 


-2.08 


± 


0.04 


-1.94 


± 


0.06 


C/2001 


RX14 (LINEAR) 


2.76 


-1.31 


± 


0.01 


-1 


39 


± 


0.00 


-1.40 


± 


0.01 


-1.36 


± 


0.00 


-1.35 


± 


0.01 


C/2002 


07 (LINEAR) (1) 


1.23 








-2 


41 


± 


0.10 


-2.49 


± 


0.28 


-2.01 


± 


0.22 


-1.06 


± 


0.34 


C/2002 


07 (LINEAR) (2) 


2.94 


-1.37 


± 


0.05 


-1 


49 


± 


0.04 


-1.47 


± 


0.04 


-1.40 


± 


0.05 


-1.37 


± 


0.05 


P/1999 


VI (Catalina) 


1.87 


-1.03 


± 


0.06 


-1 


64 


± 


0.03 


-1.60 


± 


0.03 


-1.65 


± 


0.03 


-1.55 


± 


0.10 


P/2002 


T5 (LINEAR) 


1.97 








-1 


68 


± 


0.11 


-1.70 


± 


0.08 


-1.51 


± 


0.05 


-1.34 


± 


0.53 


P/2004 


Al (LONEOS) 


0.99 








-2 


41 


± 


0.09 


-2.51 


± 


0.13 


-2.44 


± 


0.14 


-1.98 


± 


0.29 


P/2006 


U5 (Christcnsen) 


2.04 


-1.87 


± 


0.48 


-1 


57 


± 


0.05 


-1.59 


± 


0.03 


-1.70 


± 


0.03 


-1.46 


± 


0.06 



Table 3 

Linear fits of the Surface Brightness Profile slopes of the active comets in all five 
SDSS bands, described in Section 3.1.2. PSF-model is the difference between the 
SDSS measured PSF and model magnitudes 
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ct 1 


QC 


A(a)fp 


A/(0)p 


SDSS radius (km) 


Radius (km) 


Pi c f c r cn c c 


+ 1 Qp /Rnrrollv ( ~\T?C\ 
| ± y JT / DUI L\lLLy L ^ " ^ J 








2 99 


18 — 44 


h l 


Our x\.(_UI1II1LI L 11 




57 zb 2 4 


98 5 zb 4 1 


5 31 




[91 


Aft~P / Wirt Iran 




9 6 zb 2 3 


11 4 zb 2 8 


1 17 


58 


[91 




2 


26 zb 2.3 


28.3 zb 2.5 


3.56 


3.38 


[31 
L^l 


50P /Arend 


1 


12.5 zb 2.5 


18.5 ± 3.6 


2.44 


0.95 


[41 


62P /Tsuchinshan 


1 


^ 7 5 -j- 3 i 


22.7 zb 4.1 


1.09 




[21 


64P /Swift-Gchrcls 


1 


10.3 zb 1.2 


14.8 zb 1.8 


1.9 


1.83 


[21 


U 0± 1 Vj LI 11 11 \i- ) 




212 dz 22.9 


iJU ZtZ ^I.U 


8.93 


4.59 


[91 
l z J 


65P /Gunn (2) 




icq ^ 4- ■3(1 i 


&£*i2,'± l ZIZ OU.O 


8.3 


4.59 


[91 






67.2 zb 2.9 


118.4 zb 5.1 


4.84 


2.39 


[51 
l°J 


69P /Taylor 




21 zb 4.6 


31.1 zb 6.8 


1.63 


2.1 


[21 


70P /Kojima 


^ 


29.1 zb 3.5 


34.5 zb 4.2 


2.67 


2.18 


[61 
l D J 


fll3P/Spitaler (JFC) 








1.49 


1 .15 


[21 


129P /Shocmakcr-Levy 


1 


21.3 zb 2.7 


29.5 zb 3.7 


2.87 


1.66 


[21 




2 


4.4 zb 0.3 


6.8 zb 0.4 


1 .14 






1 fiSP /Tfrvwrsil T TNTT A "R ( 1 1 

1 Jor / IVUWdl- lj UN f_j r\ XV I -L J 


2 


78 3 zb 9 2 


8 ^ Q -1- 1 < 1 1 


5 42 






158P/Kowal-LINEAR (2) 


2 


98.6 zb 25.3 


110.4 zb 28.3 


6.2 






158P/Kowal-LINEAR (3) 


2 


43 zb 6.6 


50.7 zb 7.8 


5.46 






|174P (60558 Echeclus) 








39.25 


41.8 


[71 
I ' J 


fl76P (118401 LINEAR) 








2.36 






P /2002 E.I57 (LINEAR) 


2 


24.9 zb 3.8 


28.6 zb 4.4 


2.12 






2005 WY25 


2 


0.9 zb 


1.4 zb 0.1 


0.23 


0.16 


[81 
L°J 


C/1999 F2 (Dalcanton) 


1 




997.6 zb 179 


18.18 






C/2000 K2 (LINEAR) 


1 


421.6 zb 71.5 


496.5 zb 84.2 


12.36 






C/2000 Q.I46 (LINEAR) (1) 


1 


11.4 zb 0.5 


12.4 zb 0.5 


1.2 








2 


9 1 * 9-1-19 
_ ■ j . _ ZIZ J. . ^ 


31 2 zt 1 6 


3 14 






C/2000 SV74 (LINEAR) 


1 


553.1 zb 111.5 


713.2 zb 143.8 


14.63 






C/2000 Y2 (Skiff) 


2 


95.6 zb 10.4 


109.9 zb 11.9 


4.41 






C/2001 RX14 (LINEAR) 


1 


520.3 zb 53.3 


844.8 zb 86.5 


13.25 






C/2002 07 (LINEAR) (1) 


2 


152.5 zb 35.8 


187.7 zb 44 


8.87 






C/2002 07 (LINEAR) (2) 


1 


200.4 =b 18.6 


279.9 zb 25.9 


7 






C/2002 T5 (LINEAR) 


1 


92.3 zb 14.5 


112.2 zb 17.7 


6.6 






P/1999 VI (Catalina) 


1 


40.9 zb 3.9 


48.3 zb 4.7 


3.23 






P/2004 Al (LONEOS) 


2 


255.4 zb 20.1 


317.4 zb 25 


12.6 






P/2006 U5 (Christensen) 


1 


15.3 zb 1.2 


21 =b 1.7 


1.66 







Table 4 

Derived quantities. Afp values are those derived fro m SPSS observations; A(0) fp 
entri es are corrected for a phase angle of a = (see iDivind . Il98ll ; lAgarwal et al. 



2O10h . The SDSS radii are all upper limits e xcept for the unreso lved comets marked 



with a |. Literature radii values fro m [1| Lamv et al.1 ( 19981') . [2] Tancredi et all 



(1200611. [31 ISnodgrass et all (120061) . [41 Eamv et al l d200gh . [51 iTubiana et al.1 fa00$ ). 



[6] ISnodgrass et al.l &OoK [7] IStansberrv et al.l (^OoH Tand [8] Ijewittl (12006) 



Population 


u-g 


g-r 


r — i 


2 — 2 


S (u,g,r,i,z) 


S' (g,r,i) 


Comets 


1.57 ±0.21 


0.57 ±0.05 


0.22 ±0.07 


0.09 ±0.07 


8.4 ± 3.6 


8.3±3.5 


Trojans 


1.47 ±0.17 


0.58 ±0.07 


0.23 ± 0.5 


0.13 ±0.08 


8.2 ±2.9 


8.8±2.9 


Centaurs 


1.41 ± 0.46 


0.62±0.16 


0.26±0.17 


0.21±0.16 


8.9±7.7 


10.3±11.1 


TNO 


1.03 ± 1.21 


0.74 ± 0.40 


0.35 ±0.26 


0.18 ±0.17 


16.4 ±9.2 


17.8 ± 11.8 



Table 5 

Colors and spectral gradient of populations as measured by the SDSS. The colors of 
the comets, centaurs, and Trojans, and TNOs here are those of the b odies used to 



calcul ate the S , as discussed in section 4.1. The TNO colors are from lBecker et al 
(|2O10h . 
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Reference 


JFC CLF 


JFC CSD 


This Work 


0.49 ± 0.05 


-2.4 ± 0.2 


Tancredi et al. (2006) 


0.53 ± 0.05 


-2.7 ± 0.3 


Wcissman &: Lowrv (2006) 


0.35 ± 0.01 


-1.73 ± 0.06 


Lamv et al. (2004) 


0.38 ± 0.06 


-1.9 ± 0.3 


Lowrv et al. (2003) 


0.32 ± 0.02 


-1.6 ± 0.1 


Fernandez et al. (19991 


0.53 ± 0.05 


-2.65 ± 0.25 



Table 6 

Estimates of the cumulative luminosity function and cumulative size distribution 
power law slopes for the Jupiter Family Comets. 



Population 


Broken Min 


Single Power Law 


Broken Max 


Source 


fjupiter Family Comets 


0.19 ± 0.03 


0.49 ± 0.05 


0.73 ± 0.08 


This Work 


Near Earth Asteroids 




0.35 ±0.02 




Bottke et al. (2000) 


Asteroids (cometary orbits) 




0.51 ± 0.02 




Alvarez- Candal &; Licandro (2006) 


fMain Belt Asteroids (broken) 


0.10-0.59 


0.35-0.97 


0.37-1.04 


Parker et al. (2008) 


Hildas 




0.42 ± 0.02 




Alvarcz-Candal & Licandro (2006) 


f Trojans 




0.44 ± 0.05 




Szabd et al. (2007) 


Centaurs 




0.54 ± 0.07 




Larsen et al. (2001) 


Kuipcr Belt Objects 




0.66 ±0.06 




Truiillo et al. (2001) 



Table 7 

Estimates of the cumulative luminosity function power laws for small-body popu- 
lations. The CLF of Main Belt Asteroid are family populations that are fit both 
by broken and unbroken power laws, with the val ues cited here being the ranges 
determined for these families in Parker et al. ( 20081 ). Those populations marked by 
a f are derived from SDSS data. 



Comet 


R (AU) 


A/p (cm) 


Previous R (AU) 


Previous A/p (cm) 


Source 


46P/Wirtanen 


2.598 


9.59 ± 2.35 


1.120 


112.2 


1 


47P / Ashbrook- Jackson 


4.572 


25.97 ± 2.28 


4.030 


28.95 ± 4.1 


2 


65P/Gunn 


4.340 


211.97 ± 22.86 


4.430 


133.4 ± 4.7 


3 




4.338 


189.28 ± 30.37 


2.640 


23.4 


1 


f67P/ Churyumov-Gerasimenko 


1.836 


67.21 ± 2.9 


1.380 


208.9 


1 


69P/Taylor 


1.950 


21 ± 4.59 


4.030 


18.39 ± 1.7 


2 


129P/Shoemaker-Levy 


3.051 


21.26 ± 2.68 


4.558 


4.5 ± 0.6 


4 



Table 8 

A comparison of ou r A fp values wi t h pre vious listed values. Referen c es: [1 ] 
A'Hearn et all (Il995h. [2] lLowrv et al.1 (120031 1. [3] lLowrv fc Fitzsimmonsl (I2001I l 
[4] lLowrv &: Fitzsimmonsl ( 2005). ft A detai l ed loo k at Afp for 67P/Churyumov- 
Gerasimenko can be found in lAgarwal et ail (|2010h . 
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Fig. 1. Color-magnitude distribution of the resolved comets with 2-a error bars. The 
plotted symbols are: O u — g (black), A g — r (red), V r — i (green), and □ i — z 
(blue). Plotted above is the histogram of the color distribution. 
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Fig. 2. Color-color plots, showing the distribution of the SDSS comets. The resolved 
Jupiter Family Comets (JFC) are in black, while non-JFC resolved comets are in 
blue. The boxes represent the range of 2 standard de viations for the reso lved comet 
distribution (red), and Jupiter Trojans (green) from ISzabo et al.l (j2007l ). a pc is the 
prin ciple component c olor in the MOC defined as a pc = 0.89(g — r)+0.45(r — i) — 0.57 
(see llvezic et al.l l200lh . The two unresolved JFCs are in orange. 
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Fig. 3. Surface brightness profile of comet 67P/Churyumov-Gerasimenko in 
u, g, r, i, z, along with the weighted least squares fit for the slope (solid black line, 
with value in each panel's upper right). One count above sky would give the surface 
brightness shown by the horizontal dashed lines. The bottom right panel shows the 
local u — g, g — r,r — i and i — z colors with respect to radius. In all cases error bars 
are 1-a. The color remains consistent with respect to radius until the local surface 
brightness becomes close to sky values (most prominently seen with u — g in blue) 
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Fig. 4. The average slope of the u,g,and r band surface brightness profiles plotted 
against the PSF— model magnitudes of the active comets. We divide the sample into 
two "quality groups" (QG1 and QG2) at PSF— model = 1.5. For the more resolved 
QG1 comets the trend is fairly flat, and is characterized by generally small errors 
in the slope fits. The less resolved QG2 comets trend toward steeper slopes, a trend 
that terminates with slopes fit to point-sources, as shown by the four stars in the 
plot. The four objects near PSF— model = (star symbols) are field stars whose 
surface brightness profiles were extracted and fit in the same manner as the comets. 
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Fig. 5. Comparison between the best fit slopes of the surface brightness profiles for 
the SDSS filters. The solid lines are the 1:1 relations. Trends between the g, r and i 
band slopes are nearly 1:1, while there is more scatter in the u and z band trends. 
The plotted points carry la error bars. 
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Fig. 6. The SDSS colors of comets are compared against various parameters, the 
lower right plot has a key to the colors in the plot, as well as showing the median 
value for each color. There is no statistically significant trend in color against any 
parameter. The slopes of these trends have been fit, and are presented in Fig 7. 
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Fig. 7. The plot of mean color gradient of SDSS colors of comets against various 

parameters shows no color dependence on them. These same parameters can be seen 

in Figure 6. The abbreviations are: 

-SB: Slope of the surface brightness profile 

-delta_r: The r-band PSF-model magnitude 

-helio: The heliocentric distance in AU at the time of observation 
-phase: The phase angle at the time of observation 
-r_mag: The r-band magnitude 

-geo: The geocentric distance in AU at the time of observation 

-peri: The perihelion distance in AU 

-inc: The orbital inclination 

-e: The orbital eccentricity 

-Afp: The calculated Afp value in cm 

-radius: The calculated radius limit in km 
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Fig. 8. Color-color diagrams of the various small body populations discussed in 4.1. 
The active comets are in black, the Trojans are in red, the Centaurs are in orange, 
and the TNO's are in green. The error bars represent one standard deviation. 
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Fig. 9. Histograms of S for the various small body populations discussed in 4.1. 
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Fig. 10. Left: the Cumulative Luminosity Function for the Jupiter Family Comets 
in this paper, with power law fits. The whole population may be fit with a gradiant 
of 0.49 ± 0.05. Fitting with a broken power law at H ~ 14.5 gives gradients of 
0.73 ± 0.08 and 0.19 ± 0.03 for the two fits, suggesting the faint population may be 
shallower than indicated by the single fit. Right: the Size Distribution Function 
for the same sample, assuming an albedo of 0.04. There is a break between the two 
measured slopes at ~ 4km. The gradients are —3.1 ± 0.3 and —1.0 ± 0.1 for the 
broken fit. Fitting the whole range gives a fit of— 2.4±0.2. 
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Fig. 11. Calculated Afp values are compared on the left to the semi- major axis (a) 
and on the right to orbital inclination of the comets. Jupiter Family Comets are in 
black. The non-JFC population shows statistically larger Afp values. Note: LPCs 
plotted at 10 4 AU are comets on unbound orbits, indicated by arrows. 
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